Core Transitions in the Breakup of Exotic Nuclei 
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An interesting physical process has been unveiled: dynamical core excitation during a breakup 
reaction of loosely bound core + N systems. These reactions are typically used to extract spec- 
troscopic information and/or astrophysical information. A new method, the extended Continuum 
Discretized Coupled Channel (XCDCC) method, was developed to incorporate, in a consistent way 
and to all orders, core excitation in the bound and scattering states of the projectile, as well as 
dynamical excitation of the core as it interacts with the target. The model predicts cross sections 
to specific states of the core. It is applied to the breakup of 11 Be on 9 Be at 60 MeV/u, and the 
calculated cross sections are in improved agreement with the data. The distribution of the cross 
section amongst the various core states is shown to depend on the reaction model used, and not 
simply on the ground state spectroscopic factors. 

PACS numbers: 24.10.Eq, 25.60.Gc, 25.70.De, 27.20.+n 



In order to study nuclei at the limit of stability one 
needs reliable nuclear reactions models that incorporate 
the relevant structure degrees of freedom in a consistent 
manner, in particular the continuum. Theories of nuclear 
reactions have been repeatedly challenged with the new 
avenue of experimental work now possible at Radioactive 
Beam Facilities. Amongst the various reactions, breakup 
occupies a prominent role. With breakup reactions, one 
tries to extract spectroscopic information Q,J3 or capture 
reaction rates of Astrophysical relevance 0> 0, 0- I n 
either case, the structure information obtained is model 
dependent and assumes a single particle description of 
the projectile as a valence nucleon attached to the ground 
state of the core. Such a simplification may be the source 
of lingering disagreements [6j . 

When charged particle detectors are coupled with 
gamma arrays, the states of the core can be disentan- 
gled. For the variety of knockout measurements now 
available 0, the data are found to contain both elastic 
breakup (diffraction) and transfer (stripping) contribu- 
tions, which are typically calculated within an eikonal 
spectator core model 8] . Only the nuclear reaction cross 
sections to particular states are computed and other ef- 
fects need to be added a posteriori, and incoherently (i.e. 
0). One uses the single particle contributions weighted 
by the composition of the initial state predicted by shell 
model Q and neglects interference between the various 
projectile components and dynamical processes with the 
target. However, one should expect that, if the energy 
necessary to excite the core is small, there will be dy- 
namical core excitation/de-excitation during the reaction 
with the target. When a loosely bound projectile cannot 
be described as a single particle state, core excited com- 
ponents are mixed in the projectile scattering states as 
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well as the ground state. These can interfere during the 
reaction and modify the single particle picture. 

Exotic systems of the type core + N where core degrees 
of freedom may play a relevant role include 8 Li, 8 B, n Be, 
17 C, 17 F and 19 C. In addition, one can expect to find in 
their spectra states built on excited states of the core. 
This is the case of a resonance in 11 Be which is visible at 
approximately 3.4 MeV excitation in the spectrum of 0, 
and such states cannot be understood within the single 
particle description 0- This calls for a formulation of 
breakup where core excitation is consistently included. 

As a first example, we look at 9 Be( 11 Be, 10 Be)X at 60 
MeV/u Q, a reaction already studied in detail. The 
eikonal model predictions for the total cross sections to 
particular states are too low: Oth = 165 mb to be com- 
pared with er cxp = 203(31) mb for the ground state of 
10 Be, and <7 t h = 9 mb to be compared with er oxp = 16(4) 
mb for the 2 + excited state. Estimates of a Coulomb 
contribution and inelastic core excitation presented in Q 
provide a possible explanation for the apparent under- 
prediction of theory, but the problem has been await- 
ing a consistent formulation. Furthermore, Continuum 
Discretized Coupled Channel calculations show that the 
eikonal approximation does not have the desired level of 
accuracy at this energy [l(j. We propose a model where 
all these corrections and effects are included in a consis- 
tent manner. 

A recent preliminary study [Tlj generalizes the spec- 
tator core model of 8] to include core excitation. Al- 
though the work is performed for nuclear only and within 
a straight line approximation, it represents an impor- 
tant improvement over previous efforts because a core 
degree of freedom is introduced consistently in the pro- 
jectile and the core-target interaction, thus allowing for 
dynamical core excitation/de-excitation throughout the 
reaction. The initial results in pdj show an increase of 
the breakup component of the total cross sections, but 
no effect on the stripping part. This is an extremely use- 
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ful result, that the stripping component seems to be less 
affected by the various mechanisms discussed above. We 
will focus here on elastic breakup (diffraction dissocia- 
tion) only. 

As mentioned above, typically breakup models assume 
a single particle description for the projectile. Only 
recently are improvements on this approximation be- 
ing considered. One impressive improvement consists 
of describing the projectile as a three body system 
|12| adequate for nuclei of Borromean nature such as 
6 He. We pursue an alternative improvement, which is 
to describe the projectile as a multi-component system 
based on several core states. In this work, we present 
the eXtended Continuum Discretized Coupled Channel 
(XCDCC) method, to take into account explicitly core 
excitation in the breakup reaction of loosely bound sys- 
tems. 

We consider the breakup of a projectile (P), composed 
of a core (c) plus a valence particle (v), on a target 
(T). The breakup process is described using a three body 
Hamiltonian, with core degrees of freedom denoted by £: 
H 3h = T R + H proj (r, + U vt {r, R) + U ct {r, R, £). The 
coordinates are illustrated in Fig. ^ an d, as Jacobi coor- 
dinates, can be used for the full three-body wavefunction 



tions which, since the state of the core can change, are 
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where L is the projectile-target orbital angular momen- 
tum, J t the total spin of the target, Jt the total spin of 
the three-body system, and a = {L, Jp, J, J t , i, n}. 

The projectile states $j p (t*,£) can be either a bound 
or a scattering state, with several components. They are 
obtained as coupled channels eigenstates of the projectile 
Hamiltonian H pTOJ =T r + V vc {r,t) + h COTe (£): 



<:n(r)[[Yi(r)® X s] j ®MOW , (2) 



expanded in terms of core eigenstates at energies Ej. 

A continuum projectile state is characterised by 
{ Jp, i, 7i }, where i refers to the asymptotic energy, and n 
denotes the channel with a plane wave component. It is 
composed of projectile radial wavefunctions f a :n{T;k an ) 
for each a = {Isj, 1} that are solutions of coupled equa- 
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where the coupling potentials are matrix elements of 
V vc (r,£). From these equations we obtain the S-matrix 

Q 

Given the importance of continuum-continuum cou- 
plings in the CDCC formulation of breakup for halo 
nuclei [lflj . we need to transform the projectile scat- 
tering states into square integrable functions, otherwise 
continuum-continuum couplings would diverge. We de- 
fine a coupled channel bin as 



where k i 

h 2 — k 2 — 1 



I dfc e " fatnij") k a n) 
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(4) 

the core- valence relative momentum, 
■[e a — £ n ]- From S a - n we obtain S n (k), the 



FIG. 1: Coordinates for the three body breakup reaction. 



diagonal valence-core phase shift of channel a — n. Cou- 
pled channel bins defined in this way are complex. 

Substituting the three body wavefunction 
Eq. (JJJ into the three body Schrodinger equation 
H 3h ^j T = EV Jt , one obtains an equation set similar 
to the standard CDCC equations of 14], with coupling 
potentials U J J a ,{R) = (a\U ct (r,R,Z) + U vt (r,R)\a') 
containing both Coulomb and nuclear interactions 
between the projectile fragments and the target. The 
techniques for solving this equation arc the same as 
in However there are essentially two differences 

in these extended CDCC equations: first they span 
a larger number of projectile coupled components, 
corresponding to core excitation, and second the inter- 
action core-target depends on core degrees of freedom 
£. In our calculations, the couplings U a 7 a ,(R) are 
further expanded in multipoles and non-trivial algebra 
is necessary to simplify the problem, but we leave the 
details on the evaluations of these matrix elements 
for a longer publication |15| . A new generation of the 
code FRESCO 16] was developed to incorporate these 
aspects, namely coupled channel bins and core excitation 
CDCC couplings. The code was further optimized and 
parallelized so that realistic calculations could become 
feasible. 

The first application of XCDCC is to the breakup com- 
ponents of 9 Be( n Be, 10 Be)X at 60 MeV/u 0. As the spin 
s of the neutron has an insignificant effect, we set it to 
s = 0. We include only the first excited state of the core 
as m For our case the ground state of n Be contains 
two components: an s-wave neutron coupled to lo Be(0 + ) 
and a d-wave neutron coupled to 10 Be(2 + ). The model 
for n- 10 Be is based on ^8|, but the depths of the Woods- 
Saxon interactions need to be adjusted to give a positive 
parity bound state energy at E + = —0.5 MeV and a 
negative parity level at E- = —0.18 MeV, in the s = 
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model 


interaction 


V(+) 


v(-) 


R 


a 


Pi 


CC 


Be{Cr , 2^}+n 


55.25 


47.00 


2.483 


0.65 


0.67 


SP 


10 Be{0+}+n 


55.50 


30.48 


2.736 


0.67 


0.0 


SP 


10 Be{2+}+n 


75.07 


39.95 


2.736 


0.67 


0.0 


model 


interaction 


V 


W 


r 


a 


ft 


CC 


iu Be- y Be 


134.1 


68.15 


0.75 


0.8 


0.67 



TABLE I: Potential parameters for 10 Be+n with and without 
deformation of the core. Also given are the parameters for the 
core-target optical potential. 
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approximation. The needed interaction in the coupled 
channel model (CC) of n Be is listed in the first row of 
table n It produces a ground state of 11 Be where the 
neutron is 88.3% in the s-wave and 11.7% in the d-wave. 
Note that these probabilities can be related to the spec- 
troscopic factors in the shell model which are normalized 
to particle number [l7|. For comparison we also want 
to calculate the breakup cross section using the previous 
CDCC approach. In this calculation, each single particle 
contribution is multiplied by the corresponding projec- 
tile's probability and all components are added incoher- 
ently. Therefore we label this calculations by SPIS for 
single particle incoherent sum). SPIS neglects core ex- 
citation during breakup, and uses the single particle po- 
tentials built on the g.s. of the core and its first excited 
state as listed in the second and third rows of table Q] 
These potentials are used to generate the bound states 
and the whole continuum (these structure calculations 
are referred to as SP for single particle). 

The neutron-target optical potential parameters are 
taken from [l^. As the target 9 Be is a s pec tator, its 
spin J t is neglected, as in Refs [J H IIEI Oil ■ For the 
10 Be- 9 Be we take the potentials from _20j that reproduce 
unpublished data for the elastic scattering of 10 Be on 12 C 
at 59 MeV/u. This potential is directly used in the SP 
calculations. For the XCDCC calculations, we deformed 
both the nuclear and the Coulomb parts using the defor- 
mation lengths consistent with our structure model for 
n Be |2l| and refit the potential to reproduce the same 
elastic distribution. The resulting potential is given in 
the last row of table [I] 

We briefly describe the model space for solving the 
XCDCC equations. We take partial waves for Jp < 4 
organized in 178 bins as schematically shown in Fig. [21 
The evaluation of the couplings U^ T a , (R) involve an inte- 
gration in r which is performed up to r^ m = 70 fm, and a 
multipole expansion which we truncate at A max = 4. The 
projectile-target relative angular momentum is taken up 
to L max = 3000 and the corresponding distorted waves 
are matched to Coulomb functions at i? aS ym = 500 fm. 
The XCDCC calculations were performed on a SGI Altix 
3700. 

The results are shown in table |nj Using a single par- 
ticle model for 11 Be, and introducing ground state oc- 
cupation probabilities consistent with the ones produced 
by the CC model, namely 88.3 % s-wave and 11.7 % d- 
wave, the predicted cross sections are 109 mb to the g.s. 



FIG. 2: n Be continuum model space. The number of bins and 
the energy range are given for each outgoing channel (l,I n ) 
for each spin parity combination of the projectile (Jp). 



Model 


°"o+ 


cr 2+ 


G 


Eikonal 

SPIS 
XCDCC 


105 mb 
109 mb 
109 mb 


3.4 mb 
1 mb 
8 mb 


108 mb 
110 mb 
117 mb 



TABLE II: n Be breakup cross sections for 10 Be{0+,2+} + 
n. 

of the core, and 1 mb to the first excited state (model 
SPIS). In comparison, when core excitation is included 
consistently, the XCDCC calculations predict exactly the 
same cross section to the g.s. of the core, but a large in- 
crease (a factor of 8) in the cross section to the 10 Be 2 + 
state. In both SPIS and XCDCC breakup, the proba- 
bilities for seeing the 2+ state, Pi = o^/ctot, are much 
lower than in the ground state. This reflects the fact that 
the partial cross section decreases rapidly with increasing 
single-particle Q-value for breakup. 

The occupation probability (or for that matter the 
spectroscopic factor) is not an observable, therefore it 
is not good practice to compare theory and experiment 
at this level. Rather one should include an appropriate 
reaction model and compare cross sections. If we take 
the eikonal prediction for nuclear breakup directly from 
and include occupation probabilities consistent with 
our 11 Be coupled channel model, we obtain cross sections 
lower than the XCDCC predictions (see first row of table 
In}) . The differences between the theoretical cross sections 
predicted in the Eikonal model and the data Q were at- 
tributed to the Coulomb breakup and inelastic excitation 
of the core. In XCDCC one can turn off the various cou- 
plings to understand their relative importance. We find 
that even when core excitation couplings are not included 
in the core-target interaction, the cross section to the 2 + 
is increased over the SPIS prediction, due to construc- 
tive interference between the projectile's components in 
the continuum. There are also nuclear- Coulomb interfer- 
ence effects. Consequently, inelastic and Coulomb con- 
tributions should not be added incoherently. These are 
automatically contained in the XCDCC predictions. 

For a meaningful comparison with the data one needs 
to construct the stripping component, since the neu- 
tron was not detected in the measurement 0. However 
XCDCC produces the elastic breakup component only. 
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109 mb 
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200 


203(31) 
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8 mb 


6 mb 


14 
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TABLE III: Comparison of calculated cross sections, sum of 
XCDCC breakup and Eikonal stripping, with the data |l| for 
9 Be( n Be, 10 Be)X at 60 MeV/u. 

In principle we would like to produce the stripping in the 
same framework as the breakup but this is at present not 
possible. From the work of Batham et al. we learned 
that the stripping is hardly affected by core excitation, 
so we take the stripping contribution calculated in the 
eikonal approximation |22| , based on the same optical po- 
tentials as the XCDCC calculation here presented. The 
results are summarized in table ITTT1 and immediately one 
can see that the theoretical cross sections agree perfectly 
with the data. 

Since the core-target optical potential was scaled from 
that obtained from elastic scattering of 10 Be on a simi- 
lar target ( 12 C), we test the sensitivity on the choice of 
this potential. For this purpose, we use a microscopically 
based potential, by folding the NN interaction over the 
density of the core, to reproduce the same elastic and 
inelastic cross sections for 10 Be as the potential of table 
EJ The XCDCC breakup cross sections to each core state 
are not affected, and neither are the corresponding strip- 
ping cross sections. This demonstrates that details of the 
core-target optical potential are not important, as long 
as similar observables for the core-target interaction are 
obtained. 

Finally, it is necessary to point out that within our 
model space, there is no population of the 10 Be 1~ and 
2~ states, in this reaction. These states are seen in the 
experiment Q, but we expect that these states, which 
result from neutron excitation from the core, will not 
interfere with the results here presented. At present we 
do not have a 10 Be structure model that incorporates all 
2+,l",2~ states in a simple way, as these correspond to 
breaking the core. 

In conclusion, we find that the amount of core exci- 
tation is modified in the breakup reaction. This process 
occurs both through constructive interference of various 
components in the projectile and through the interac- 



tion of the core with the target. The extended Contin- 
uum Discretized Coupled Channel method was developed 
specifically to handle the problem of core excitation in the 
breakup of loosely bound projectiles. Effects of core ex- 
citation in the projectile bound and scattering states are 
explicitly taken into account. Through the interaction 
with the target, the core can excite or de-excite during 
the reaction. We have applied XCDCC to the breakup 
of n Be on 9 Be at 60 MeV/u, where the final 10 Be state 
is identified. Theoretical predictions within a truncated 
model space agree with the data for the cross sections to 
the first two individual 10 Be states. 

As compared to the preliminary calculations of 
Batham et al. [HI, XCDCC provides more detailed cross 
sections, namely partial cross sections to each core state. 
As a consequence we now understand that the increase 
of the total cross section when including core excitation 
comes mainly from an increase of the core excited com- 
ponent. Thus one can think of it as production of core 
excitation. This process was not understood before. 

Despite the computational challenge, XCDCC consists 
on a significant improvement to previous theories. Given 
the possibility of producing core excitation, previous 
spectroscopic analyses and extractions of astrophysical S- 
factors need to be revisited. Other reactions that can be 
usefully studied with XCDCC include 11 Be(p,p') 11 Be*, 
12 C( 11 Be, 10 Be+n), the various modes of 8 B breakup, and 
work along these lines is underway. So far, the method is 
limited to the inclusion of states of the core that can be 
modelled as collective excitations, but it could easily be 
adapted to including a better description of the core, as 
long as a complete set of core + N (bound and scattering 
states) could be obtained. 

We thank Jeff Tostevin for providing the stripping 
cross sections for comparison with the data and for use- 
ful discussions. We also thank Jim Al-Khalili for pro- 
viding the data for 10 Be. This work is supported by 
NSCL, Michigan State University, the U.K. EPSRC by 
grant GR/T28577, and by the National Science Founda- 
tion through grant PHY-0456656. The calculations pre- 
sented were performed at the HPCC at Michigan State 
University. 
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